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Abstract 

A three-dimensional time-dependent quantum dynam- 
ics approach using a recently developed ab initio poten- 
tial energy surface is applied to study ro- vibrational exci- 
tation in N+N 2 exchange scattering for collision energies 
in the range 2.1- 3.2 eV. State-to-state integral exchange 
cross sections are examined to determine the distribution 
of excited rotational states of N 2 . The results demon- 
strate that highly-excited rotational states are produced 
by exchange scattering and furthermore, that the max- 
imum value of A j increases rapidly with increasing col- 
lision energies. Integral exchange cross sections and ex- 
change rate constants for excitation to the lower ( v — 
0-3) vibrational energy levels are presented as a function 
of the collision energy. Excited-vibrational-state distri- 
butions for temperatures at 2,000 K and 10,000 K are 
included. 


Nomenclature 

a 0 = Bohr radius, 5.291772 xlCT 9 cm 

on — coefficient in Eq. (3) 

D = product vibrational distribution 

E = Total scattering energy 

/ = subscript denoting final state 
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h = Plank constant 

i — subscript denoting initial state 

j = N 2 rotational quantum number 

J = total angular momentum 

k = reaction rate constant, cm 3 /sec 

P = state-resolved reaction probability 

r = N-N separation distance, ao 

R = distance between N and N 2 midpoint, ao 

t = propagation time 

T — temperature, K 

v = N 2 vibrational quantum number 

A j = change in value of j 

k = Boltzmann constant 

fi = reduced mass 

a — integral cross section, Uq 

6 = angle between r and R 

0 = angle between r a and rt, 

ft = projection of J along z axis 

T = time-dependent wave function 

= time-independent wave function 


Introduction 

In a recent shock tube experiment using nitrogen gas, 
Fujita et al. 1 measured the rotational and vibrational 
temperatures using the N 2 2nd postive and first 
negative bands and found significant nonequilibrium be- 
tween the measured rotational and translational temper- 
atures. This result is in apparent violation of the two- 
temperature model of Park 2 , one of the standard models 
widely used in the analysis of nonequilibrium aerother- 
modynamics. The observed deviation led to two studies 
of nitrogen rotational excitations behind a strong shock 
wave. Fujita and Abe 3 carried out quasi-classical tra- 
jectory (QCT) calculations of state-to-state rotational 
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FIG. 1: Contour plot of the WSHDSP ab initio PES of 
the N+N 2 interaction region in terms of the two N-N bond 
lengths r a and r*>; © is fixed at 120°. 
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FIG. 2: Contour plot of the LEPS PES for the N+N 2 inter- 
action region of Fig 1 in terms of the two N-N bond lengths 
r a and rb\ 0 is fixed at 180°. 


excitation rates in N2T-N2 collisions using the empiri- 
cal N2-N2 intermolecular potential energy of Billing and 
Fisher 4 . Their calculation covered transitions with very 
large (well over 100) A j, which is the change in the N2 
rotational state quantum number j .Their QCT rate dif- 
fered significantly from the empirical Rahn-Palmer 5 rate 
formula above 3000 K. Simultaneously, Park 6 extended 
the Rahn-Palmer formula to large A j and high tempera- 
ture. His results appeared to agree reasonably well with 
the shock tube data of Sharma and Gillespie 7 and Ref. 
1 . 

Both approaches, the calculation of Ref. 3 and the 
Rahn-Palmer extension of Ref. 6 have substantial short- 
comings that make their predictions unreliable. St allcop 
and Partridge 8 have constructed a potential energy sur- 
face (PES) for N2--N2 interactions that is based on ac- 
curate short-range ab initio results from high-level quan- 
tum structure calculations; the van der Waals interac- 
tion and the long-range behavior of this PES have been 
validated 9 by very accurate results from a bond-function 
calculation. Knowing the anisotropy of the actual N2-N2 
interaction energy assists our assesment of the reliability 
of the two approches. 

The N2-N2 interaction energy used in calculation of 
Ref. 3 was developed in Ref. 4 to study V-V and V-T 


energy transfer rates. Rotational motion was not quan- 
tized in Ref. 4, but was considered to be part of the 
translational motion. The angular dependence of the re- 
pulsive interaction energy was represented by the factors: 
cosh(ei) with c* = 0.5aricos^i, where r* is the N-N sep- 
aration distance of the fth N 2 molecule and 7 \ the angle 
between r* and R, which is the separation distance of the 
center of masses of the two N2 molecules. Also, a is a 
parameter that specifies the the range of the repulsive 
potential energy. The use of the cosh(ei) factors is un- 
usual. Because the repulsive interaction energy of Ref. 
4 does not fully reflect the the symmetry of the N2-N2 
interaction, it can lead to very high A j transitions that 
are physically unrealistic. It has been shown that the 
angular variation of the repulsive potential energy of the 
accurate interaction energy can be represented by a rel- 
atively simple form 8 , but it is very different from that 
of Ref. 4. Because rotational motion was not explicitly 
accounted for in the calculation of Billing and Fisher, the 
above factors did not create readily-discernable problems 
in their results 4 . On the other hand, rotational motion 
was quantized in the study of Ref. 3. The cosh(ci) factors 
were the source of very large A j transitions, which would 
not have been found if a more realistic N2-N2 interaction 
energy had been used. 




FIG. 3: State-to-state integral cross sections for exchange 
scattering, with rotational state excitation, N + N 2 (u = 0,j — 
0) — ► N 2 (*/ = 0,/) 4- N, vs. j for total scattering energy at 
2.5 eV. The symbols o and o represent the results obtained by 
neglecting nuclear spin and using Bose-Einstein statistics, re- 
spectively. The discrete data points are joined by straight-line 
curves for clarity. 

The extension of the Rahn-Palmer formula used in Ref. 
6 to predict rotational excitation at large A j and high 
temperatures probably also produces unphysical results. 
The Rahn-Palmer formula is fitted to experimental data 
only up to 3000 K; the largest observed A j was 14. Ex- 
trapolation to very large A j (such as more than 100) 
is beyond the range of validity of this formula. Based 
on the actual N 2 -N 2 interaction energy described above, 
we conclude that such highly anisotropic interaction is 
probably unphysical. 

While it is likely that the large A j transitions indicated 
in the shock tube data of Ref. 1 do not come from N 2 - 
N 2 collisions, the N-N 2 exchange mechanism can lead to 
such transitions at the high temperature regime of the ex- 
periment. Exchange reactions are important phenomema 
for real-gas analysis. For example, like ion-atom diffusion 
processes are dominated by resonance charge exchange. 

In the atom exchange reaction 

N+N2(vJ)^N 2 (v'J') + N, ( 1 ) 

the N-N bond in the original N 2 molecule is broken 
and a new N-N bond is formed. Due to the nature of 
the N+N 2 interaction energy as described below, the 
new N 2 molecule may have a highly-excited rotational 
state (compared to the rotational state of the original 
molecule). 



FIG. 4: Similar to Fig. 3 except that the total scattering 
energy is at 3.0 eV. 


Among more recent model calculations, Esposito et 
al. 10 ’ 11 have calculated N+N 2 scattering using a quasi- 
classical method and a semiempirical London- Eyring- 
Polany-Sato (LEPS) interaction energy. All but one of 
parameters of the LEPS model have been obtained 12 
from the potential energy well of a ground-state N 2 
molecule or an extrapolation using a Morse-potential- 
energy fit; the remaining parameter specifies the height 
of a transition state barrier and has been estimated from 
the results of limited measurement. 

A theoretical determination of exchange scattering 
must be based on a quantum mechanical description. For 
example, the N-N + resonance charge exchange 13 would 
be underestimated by almost an order of magnitude using 
classical mechanics. The first ab initio potential energy 
surface for the N 3 system (WSHDSP PES) that includes 
variation in the interatomic separation distances of the 
N 2 molecules and the first quantum dynamics study of 
the cumulative reaction probabilities and chemical reac- 
tion rates of the N-{-N 2 exchange reaction have recently 
been reported by our group 14 . 

The WSHDSP 14 PES for N 3 is constructed from the re- 
sults of high-level quantum structure calculations. These 
calcuations axe primarily based on the same methods 
that were used to determine the ab initio data for the 
construction of the very accurate N+ N 2 rigid-rotor PES 
of Stallcop et al. 15 ; long-range dispersion force data and 
nearly 4000 ab initio data points, are applied to construct 
our present N 3 PES. The reaction region of the PES con- 





FIG. 5: State to state integral cross sections for exchange 
scattering with with rotational state excitation, N + N 2 (v = 
0, j — 0) — ► N 2 (v' — 0 ,/) + N, as a function of the kinetic 
energy. 


tains two transition state barriers, which are equivalent 
with respect to interchange of two nitrogen atoms, that 
are separated by a shallow well. The ab initio PES 14 
and the results from a semiempirical model 12 are repre- 
sented by contour plots in Figs. 1 and 2, respectively. 
Note the large difference in the behaviors of the actual 
and semiempirical PES in the reaction region. 

The reaction geometries of our PES is specified by: 
r a (ri) = 2.23 bohr, n(r a ) = 2.80 bohr, and 0 — 119° for 
the transition states and r a — = 2.40 bohr and 0 = 

120° for the well minimum; r a and are the two N-N 
bond lenths and 0 is the angle between r a and r&. The 
transition state is located at 2.05 eV above the energy of 
the N-hN 2 asymptote, and the well minimum is at 1.90 
eV. The shallow well of the double hump in the reaction 
region has been dubbed 16 * 17 “Eyring’s lake” in studies 
of other interactions; in the discussion below, we use the 
term 'Lake Eyring’ to identify certain features associated 
with reaction region of the N3 PES. 

In addition to the attributes noted above, a quantum 
mechanical description of the scattering can account for 
inelastic scattering reasonances and effects due to nuclear 
spin symmetry. We have found 18,19 that N+N 2 exchange 
scattering exhibits profuse Feshbach resonance structures 
and furthermore, that these resonances arise from N3 vi- 
brational energy levels lying above the shallow well of the 
transition-state region of the PES. 



FIG. 6: Integral cross sections for exchange scattering with 
N + N 2 (v = 0, j = 0) — > Naiv') -h N as a function of kinetic 
energy; v' ~ 0-3. 

Theoretical Method 

The 3D quantum dynamics of the state-to-state reac- 
tion calculation is carried out using a time-dependent 
wave packet method 20,21 ; this method is summarized 
briefly below. The split-operator 22 propagation scheme 
is used for wave packet propagation. Also, a two-stage 
propagation approach is employed to calculate the state- 
resolved reaction probability. First, an initial wave 
packet for a given total angular momentum J, located in 
the reactant asymptotic region and covering the transla- 
tional energy from 2.0 eV to 3.2 eV, is propagated in the 
reactant Jacobi coordinates on the ab initio PES, until 
it enters into the transition state region after about 2000 
atomic units of time. At this stage, its position represen- 
tation is transformed from the reactant into the product 
Jacobi coordinates and the propagation is continued. 

After the wave packet enters into the product asymp- 
totic region, the energy-dependent state-resolved reac- 
tion probability Pfi(E) for energy E is obtained from 
the flux passing through a dividing surface in the prod- 
uct asymptotic region 

Pfi(E) = (2) 

where p is the reduced mass of the system. The quan- 
tities i and / represent the set of quantum numbers 
v,j, J, O and v\j\ J, fT for the initial and final states, 
respectively. Also, R is the translation Jacobi coordinate 
that connects an N atom to the center mass of the N 2 



molecule; unprimed or primed indicates that the coordi- 
nate system is located in the reactant or product channel, 
respectively. The quantity R'j is the position of the di- 
viding surface in the product asymptotic region where 
the reaction is complete and the flux is invariant with 
respect to moving Rj further away from the interaction 
region. The time-independent scattering wavefunction 
R r ) is obtained by performing a Fourier transform 
of the time-dependent wave equation t.e., 

1 /»+ OO 

^ iiE ' R ' ) = ^(E)J „ e <Tt </.-RW«)>*- (3) 

where the coefficient ai(E) is given by the overlap 23 be- 
tween the initial wave packet 4q(0) and the normalized 
asymptotic scattering function. 

Since a two-stage propagation scheme is employed to 
calculate the state-to-state reaction probability, there are 
two sets of basis functions for this calculation. One set 
is for the reactant channel, with 220 sine functions to 
represent R from 2.3 ao to 10.4 ao, 100 vibrational eigen- 
functions of N 2 to cover the vibrational coordinate (N-N 
separation distance) r from 1.3 ao to 5.5 ao, and 70 even 
Legendre functions (j ma x=l 40) or 70 odd Legendre func- 
tions (imax= 139) to cover 0, which is defined to be the 
angle between R and r. The other set is for the prod- 
uct channel, where 320 sine functions are used for R' in 
the range from 0.1 ao to 14.0 ao, 220 vibrational func- 
tions for expansions of r' from 0.1 ao to 12.0 ao, and 150 
Legendre functions for & . As described above, the wave 
packet first propagates in the reactant channel and then 
it is transformed into the product-channel coordinates. 
Finally, it is propagated for 8000 atomic units of time in 
the product channel. The state resolved flux is calculated 
at Ry— 10.0 ao. The whole calculation has been carried 
out using 32 CPU on a SGI Origin 3000 system. About 
40,000 CPU hours were required for the state-resolved 
calculations described here. 

When nuclear spin is ignored (no-spin), the integral 
cross section can be obtained by summing over all partial 
waves J 

= kfWT T) £ (2J + 1)p ^ ( £)- (4) 

where P^ j/4 _ v j(E) is obtained from the Pfi of Eq. (1) 
by summing over Q and Q'. The quantity ki in Eq. (4) 
comes from the initial-state wave vector. 



FIG. 7: Exchange rate constant for N + N 2 (u = 0 t j = 0) — ► 
N 2 (t/) + N; T is in K, k is in cm 3 /sec, and v' = 0-3. 


Since the nitrogen atoms have identical nuclei (99.7% 
of nitrogen consists of the 14 N isotope), the N-fN 2 cross 
section a is obtained from sums and averages over the 
nuclear spin states. The nuclear spin of nitrogen atom 
is equal to 1; hence, the N-fN 2 scattering obeys Bose- 
Einstein (BE) statistics. The cross sections with the ap- 
propriate nuclear spin factors become 


<?v' j' <—vj — 



if j' is even, 
if f is odd. 


( 5 ) 


Results and Discussion 
Rotational Excitation 

The determination of the state-to-state integral cross 
section requires summing over 96 partial waves ( i.e ., from 
J — 0 to Jmax = 95) to obtain convergence for kinetic 
energies in the range from 2.0 eV to 3.2 eV. State-to- 
state integral cross sections for scattering energies at 2.5 
eV and 3.0 eV are given in Figs. 3 and 4, respectively, to 
illustrate the distribution of rotational state excitation. 
Note that we plot both results, one with nuclear spin 
based on Bose-Einstein statistics and the other without 
spin. Since the nuclear spin factors of Eq. (5) introduce 
rapid oscillations in the cross sections, resonance features 
in the cross sections are more easily recognizable using 
the data without nuclear spin. In Fig. 3, strong oscil- 
lations occur at low f and at the three major reactive 
resonances centered at f — 19, 30, and 43. There is no 




significant excitation to rotational states with f larger 
than 57. Comparison of Fig. 3 with the similar results 
shown in Fig. 10 of Ref. 19 for a scattering energy at 2.4 
eV shows that the excited rotational distribution is quite 
sensitive the scattering energy. Fig. 4 shows that higher 
(up to f = 70) rotational states of N 2 are excited at a 
higher scattering energy. Furthermore, the integral cross 
sections are larger at higher scattering energies; strong 
oscillations are also found at low f and have a variation 
similar to that shown in the figures discussed above. 

Figure 5 shows a V 'j^ V j for scattering from a (v = 
0J = (J) initial ro- vibrational state to various final ro- 
vibrational states. The curves of crof<-oo for low-lying 
final ro-vibrational states, such as the ( v' — 0 J' = 0) 
and (u' = 0j' ~ 2) states, have an oscillating behavior 
arising from reactive resonances. As described above, 
these resonances have been shown 18,18 to be caused by 
metastable states of the ‘Lake Eyring’ region of the PES. 
Here they appear, primarily, as ‘bumps’ in the state-to- 
st ate integral cross sections. However, for higher final 
ro-vibrational states, e.g ., ( v ' = 0,/ = 48) and ( v f = 
0 J' = 69) states, the resonances are washed out by the 
sum over the partial waves. Fig. 5 also shows that the 
integral cross section for j = 0 — ► f = 2 has a much 
smaller magnitude than that for j = 0 — » f = 48 at lower 
kinetic energies. The variation of the large cross sections 
for producing high /, such as j = 0 — * f = 48, with 
kinetic energy verifies that the N-j-N 2 exchange reaction 
is quite efficient for producing high rotational levels of 
N 2 molecules. 

The above results for rotational excitation indicate 
that the exchange reaction of N-bN 2 collisons provides 
a mechanism to pump N 2 molecules with low rotational 
states into N 2 molecules with higher rotational states. 
Because this exchange reaction has a propensity for pro- 
ducing highly excited rotational states, it provides a 
route that may lead to nonequilibrium rotational tem- 
perature behind a strong shock. 

Vibrational Excitation 

The integral cross section for v f excitation from a 
(n — 0,j = 0) initial ro-vibrational state is obtained by 
summing over all the final rotational states ; i.e., 

& v /<— 00 = <7Vjf/<„00* (6) 

P 



FIG. 8: The product vibrational distribution for N + N 2 (v = 
0,jf = 0) — > N 2 (u') 4- N vs. v* . The discrete data for each 
value of T (at 2000 and 10,000 K) are joined by straight-line 
curves for clarity. 

The values of cr*,^ 00 are shown in Fig. 6 for v' = 0-3. 
This figure shows that significant vibrational excitation 
of N 2 is produced by exhange reactions at higher kinetic 
energies. In particular, this figure shows that at lower ki- 
netic energies, the largest integral cross section is found 
for v' = 0. However, at higher energies the order is differ- 
ent; for example, at 3.2 eV the order is: o-*^ 00 > <r 2 <_oo > 
^ 0^-00 > f n other words, at this energy, the vibra- 

tional products are larger at higher (v' ~ 2 and t/ = 1) 
vibrational levels than at the lower (v f = 0) vibrational 
level. Furthermore, the v' = 3 vibrational state does not 
have a major role in the final state distributions at this 
energy. 

Rate Constants 

The reaction rate constant fc v /^oo(T 1 ), for exchange 
scattering from the (v — 0,j = 0) initial ro-vibrational 
state to the (v\f) final ro-vibrational states for all /, is 
obtained from 

k v ,^oo(T) = ^ 3^3 J dEexp{-E/KT)Ea* v ,^ m {E). 

( 7 ) 

We present for v' = 0-3 in Fig. 7. Note that 

the rates can be represented fairly well by straight lines 
except at high T where there is a small fall-off; the inte- 
gration of Eq. (7) has been truncated at the the highest 
energy of this study. 

The product vibrational distributions are calculated 
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from the equation 

Dv'<-oo(T) = Av^oo^V^^'^ooCF)- (8) 

V " 

Figure 8 shows A/<_ 00 (T) for T — 2,000 K and 10,000 K. 
Note that the distribution of the final vibrational states 
of N 2 is more concentrated at low tf . As T increases, the 
peak of the distribution moves to a higher-lying vibra- 
tional state as expected. 

Concluding Remarks 

Our calculations demonstrate that the N+N 2 exchange 
reaction could result in a rotational distribution that is 
very much in nonequilibrium. This was not predicted by 
a previous QCT study 10,11 of N-bN 2 collisions because 
the LEPS interaction energy 12 , which was used for their 
study, does not have the Take Eyring’ feature. Further- 
more, the resonance structures in the reaction probabil- 
ity are quantal features that are not amenable to classical 
treatment. 

One might disregard like-atom exchange from atom- 
molecule collisions since there is no change in the pop- 
ulation density of the constituents of the gas; the above 
study, however, demonstrates that these reactions could 
significantly change the energy distribution among the 
constituents. Furthermore, we point out that the above 
approach, construction of a PES from ab initio data and 


subsequent quantum-dynamical calculation of the atom- 
molecule reactions, can also be applied to determine the 
exchange reaction for unlike atoms, such as 

O + N 2 -> NO + TV. (9) 

The use of model interaction energies in scattering cal- 
culations may reproduce the results of certain experi- 
ments, but we emphasize that, in general, the results of 
these calculations do not provide reliable thermodynamic 
properties, especially at high temperatures. 

Direct N+N 2 inelastic collisions, which do not involve 
any bond breakage and formation, will lead to low A j 
transitions until the kinetic energy is at or above the 
onset of exchange reaction. In this region, the resonances 
that play such an important role in exchange collisions 
also affect inelastic collisions and large A j transitions are 
expected from these inelastic collisions, as well as that 

found for the exhange reactions. 

We plan to study ro-vibrational excitations in N+N 2 

inelastic collsions. The inelastic results from direct scat- 
tering will then be incorporated with the present results 
for N+N 2 exchange collisions. The data will be used to 
similate the shock tube experiments of Refs. 1 and 7. Ro- 
tational and vibrational temperatures under high-speed 
entry conditions will also be studied. 
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